Rational Design, Synthesis and Binding Affinity Studies of Anthraquinone Derivatives Conjugated to Gonadotropin-Releasing Hormone (GnRH) Analogues towards Selective Immunosuppression of Hormone-Dependent Cancer.

Gonadotropin-releasing hormone (GnRH) is pivotal in regulating human reproduction and fertility through its specific receptors. Among these, gonadotropin-releasing hormone receptor type I (GnRHR I), which is a member of the G-protein-coupled receptor family, is expressed on the surface of both healthy and malignant cells. Its presence in cancer cells has positioned this receptor as a primary target for the development of novel anti-cancer agents. Moreover, the extensive regulatory functions of GnRH have underscored decapeptide as a prominent vehicle for targeted drug delivery, which is accomplished through the design of appropriate conjugates. On this basis, a rationally designed series of anthraquinone/mitoxantrone-GnRH conjugates (con1-con8) has been synthesized herein. Their in vitro binding affinities range from 0.06 to 3.42 nM, with six of them (con2-con7) demonstrating higher affinities for GnRH than the established drug leuprolide (0.64 nM). Among the mitoxantrone based GnRH conjugates, con3 and con7 show the highest affinities at 0.07 and 0.06 nM, respectively, while the disulfide bond present in the conjugates is found to be readily reduced by the thioredoxin (Trx) system. These findings are promising for further pharmacological evaluation of the synthesized conjugates with the prospect of performing future clinical studies.

Gonadotropin-Releasing Hormone and GnRH Receptor: Structure, Function and Drug Development.

BACKGROUND: Gonadotropin-Releasing Hormone (GnRH) is a key element in sexual maturation and regulation of the reproductive cycle in the human organism. GnRH interacts with the pituitary cells through the activation of the Gonadotropin Releasing Hormone Receptors (GnRHR). Any impairments/dysfunctions of the GnRH-GnRHR complex lead to the development of various cancer types and disorders. Furthermore, the identification of GnRHR as a potential drug target has led to the development of agonist and antagonist molecules implemented in various treatment protocols. The development of these drugs was based on the information derived from the functional studies of GnRH and GnRHR. OBJECTIVE: This review aims at shedding light on the versatile function of GnRH and GnRH receptor and offers an apprehensive summary regarding the development of different agonists, antagonists and non-peptide GnRH analogues. CONCLUSION: The information derived from these studies can enhance our understanding of the GnRH-GnRHR versatile nature and offer valuable insight into the design of new more potent molecules.

Refinement of the gonadotropin releasing hormone receptor I homology model by applying molecular dynamics.

Sexual maturation of human cells in ovaries and prostate is linked to the biochemical cascade initiated by the activation of cell receptors through the binding of Gonadotropin Releasing Hormone (GnRH). The GnRH receptors (GnRHR) are part of the rhodopsin G-protein coupled receptor (GPCR) family and consist of seven trans-membrane helical domains connected via extra- and intra-cellular segments. The GnRH-GnRHR complex has been implicated in various forms of prostate and ovarian cancer. The lack of any structural data about the GnRH receptor impedes the design of antagonists for use in cancer treatment. The aim of the study is to devise a model of GnRHR to be used further for the design of improved peptide/non-peptide GnRH analogues and, to our knowledge provide new structural information regarding the extracellular loop 2 (ECL2) that acts a regulator of ligand entry to GnRHR. The common structural characteristics, of the members of the rhodopsin family of GPCRs, have been employed for the construction of a homology model for GnRHR. Structural information from the human ß2-adrenergic receptor, as well as rhodopsins have been used in order to create a theoretical model for GnRHR. Furthermore, molecular dynamics (MD) simulations have been employed for the refinement of the model and to explore the impact of the bilayer membrane in GnRHR conformation.

In Silico Drug Design: Non-peptide Mimetics for the Immunotherapy of Multiple Sclerosis.

Advances in theoretical chemistry have led to the development of various robust computational techniques employed in drug design. Pharmacophore modeling, molecular docking, and molecular dynamics (MD) simulations have been extensively applied, separately or in combination, in the design of potent molecules. The techniques involve the identification of a potential drug target (e.g., protein) and its subsequent characterization. The next step in the process comprises the development of a map describing the interaction patterns between the target molecule and its natural substrate. Once these key features are identified, it is possible to explore the map and screen large databases of molecules to identify potential drug candidates for further refinement.Multiple sclerosis (MS) is an autoimmune disease where the immune system attacks the myelin sheath of nerve cells. The process involves the activation of encephalitogenic T cells via the formation of the trimolecular complex between the human leukocyte antigen (HLA), an immunodominant epitope of myelin proteins, and the T-cell receptor (TCR). Herein, the process for rational design and development of altered peptide ligands (APLs) and non-peptide mimetics against MS is described through the utilization of computational methods.

Biologically relevant conformational features of linear and cyclic proteolipid protein (PLP) peptide analogues obtained by high-resolution nuclear magnetic resonance and molecular dynamics.

Proteolipid protein (PLP) is one of the main proteins of myelin sheath that are destroyed during the progress of multiple sclerosis (MS). The immunodominant PLP139-151 epitope is known to induce experimental autoimmune encephalomyelitis (EAE, animal model of MS), wherein residues 144 and 147 are recognized by T cell receptor (TCR) during the formation of trimolecular complex with peptide-antigen and major histocompability complex. The conformational behavior of linear and cyclic peptide analogues of PLP, namely PLP139-151 and cyclic (139-151) (L144, R147) PLP139-151, have been studied in solution by means of nuclear magnetic resonance (NMR) methods in combination with unrestrained molecular dynamics simulations. The results indicate that the side chains of mutated amino acids in the cyclic analogue have different spatial orientation compared with the corresponding side chains of the linear analogue, which can lead to reduced affinity to TCR. NMR experiments combined with theoretical calculations pave the way for the design and synthesis of potent restricted peptides of immunodominant PLP139-151 epitope as well as non peptide mimetics that rises as an ultimate goal.

Design and Synthesis of Non-Peptide Mimetics Mapping the Immunodominant Myelin Basic Protein (MBP83-96) Epitope to Function as T-Cell Receptor Antagonists.

Encephalitogenic T cells are heavily implicated in the pathogenesis of multiple sclerosis (MS), an autoimmune demyelinating disease of the central nervous system. Their stimulation is triggered by the formation of a trimolecular complex between the human leukocyte antigen (HLA), an immunodominant myelin basic protein (MBP) epitope, and the T cell receptor (TCR). We detail herein our studies directed towards the rational design and synthesis of non-peptide mimetic molecules, based on the immunodominant MBP83-96 epitope that is recognized by the TCR in complex with HLA. We focused our attention on the inhibition of the trimolecular complex formation and consequently the inhibition of proliferation of activated T cells. A structure-based pharmacophore model was generated, in view of the interactions between the TCR and the HLA-MBP83-96 complex. As a result, new candidate molecules were designed based on lead compounds obtained through the ZINC database. Moreover, semi-empirical and density functional theory methods were applied for the prediction of the binding energy between the proposed non-peptide mimetics and the TCR. We synthesized six molecules that were further evaluated in vitro as TCR antagonists. Analogues 15 and 16 were able to inhibit to some extent the stimulation of T cells by the immunodominant MBP83-99 peptide from immunized mice. Inhibition was followed to a lesser degree by analogues 17 and 18 and then by analogue 19. These studies show that lead compounds 15 and 16 may be used for immunotherapy against MS.

Molecular dynamics at the receptor level of immunodominant myelin oligodendrocyte glycoprotein 35-55 epitope implicated in multiple sclerosis.

Multiple Sclerosis (MS) is a common autoimmune disease whereby myelin is destroyed by the immune system. The disease is triggered by the stimulation of encephalitogenic T-cells via the formation of a trimolecular complex between the Human Leukocyte Antigen (HLA), an immunodominant epitope of myelin proteins and T-cell Receptor (TCR). Myelin Oligodendrocyte Glycoprotein (MOG) is located on the external surface of myelin and has been implicated in MS induction. The immunodominant 35-55 epitope of MOG is widely used for in vivo biological evaluation and immunological studies that are related with chronic Experimental Autoimmune Encephalomyelitis (EAE, animal model of MS), inflammatory diseases and MS. In this report, Molecular Dynamics (MD) simulations were used to explore the interactions of MOG35-55 at the receptor level. A detailed mapping of the developed interactions during the creation of the trimolecular complex is reported. This is the first attempt to gain an understanding of the molecular recognition of the MOG35-55 epitope by the HLA and TCR receptors. During the formation of the trimolecular complex, the residues Arg(41) and Arg(46) of MOG35-55 have been confirmed to serve as TCR anchors while Tyr(40) interacts with HLA. The present structural findings indicate that the Arg at positions 41 and 46 is a key residue for the stimulation of the encephalitogenic T-cells.

Elucidation of the binding mechanism of renin using a wide array of computational techniques and biological assays.

We investigate the binding mechanism in renin complexes, involving three drugs (remikiren, zankiren and enalkiren) and one lead compound, which was selected after screening the ZINC database. For this purpose, we used ab initio methods (the effective fragment potential, the variational perturbation theory, the energy decomposition analysis, the atoms-in-molecules), docking, molecular dynamics, and the MM-PBSA method. A biological assay for the lead compound has been performed to validate the theoretical findings. Importantly, binding free energy calculations for the three drug complexes are within 3 kcal/mol of the experimental values, thus further justifying our computational protocol, which has been validated through previous studies on 11 drug-protein systems. The main elements of the discovered mechanism are: (i) minor changes are induced to renin upon drug binding, (ii) the three drugs form an extensive network of hydrogen bonds with renin, whilst the lead compound presented diminished interactions, (iii) ligand binding in all complexes is driven by favorable van der Waals interactions and the nonpolar contribution to solvation, while the lead compound is associated with diminished van der Waals interactions compared to the drug-bound forms of renin, and (iv) the environment (H2O/Na(+)) has a small effect on the renin-remikiren interaction.

Systematic molecular dynamics, MM-PBSA, and ab initio approaches to the saquinavir resistance mechanism in HIV-1 PR due to 11 double and multiple mutations.

Mutations in the human immunodeficiency virus (HIV) enable virus replication even when appropriate antiretroviral therapy is followed, thus leading to the emergence of drug resistance. In a previous work, we systematically examined seven single mutations that are associated with saquinavir (SQV) resistance in HIV-1 protease (Tzoupis, H.; Leonis, G.; Mavromoustakos, T.; Papadopoulos, M. G. J. Chem. Theory Comput. 2013, 9, 1754-1764). Herein, we extend our analysis, which includes seven double (G48V-V82A, L10I-G48V, G48V-L90M, I84V-L90M, L10I-V82A, L10I-L63P, A71V-G73S) and four multiple (L10I-L63P-A71V, L10I-G48V-V82A, G73S-I84V-L90M, L10I-L63P-A71V-G73S-I84V-L90M) SQV-HIV-1 PR mutant complexes, in an attempt to generalize our findings and formulate the main elements of the SQV resistance mechanism in the protease. On the basis of molecular dynamics (MD), molecular mechanics Poisson-Boltzmann surface area (MM-PBSA), and ab initio computational approaches, we identified specific features that constitute the HIV-1 PR mechanism of resistance at the molecular level: the low flexibility of SQV in the binding cavity and the preservation of hydrogen bonding (HB) and van der Waals interactions between SQV and several active-site (Gly27/27', Asp29/29'/30/30', especially Asp25/25') and flap (Ile50/50', Gly48/48') residues of the protease contribute significantly to efficient binding. The total enthalpy loss in all mutants is mostly due to the loss in enthalpy of the active-site region. Furthermore, it was observed that mutation accumulation may induce stabilization to SQV and to the flaps through enhanced HB interactions that lead to improved inhibition (e.g., accumulation of mutations in complexes containing L10I, G48V, L63P, I84V, or L90M single mutations). It was also concluded that permanent flap closure is obtained independently of mutations and SQV binding is mostly driven by van der Waals, nonpolar, and exchange-energy contributions. Importantly, it was indicated that the optimal positioning of SQV and the structure of the binding cavity are tightly coupled, since small changes in geometry may affect the binding energy greatly. The results of our theoretical approaches are in agreement with experimental evidence and provide a reliable description of SQV resistance in HIV-1 PR.

A Comparative Molecular Dynamics, MM-PBSA and Thermodynamic Integration Study of Saquinavir Complexes with Wild-Type HIV-1 PR and L10I, G48V, L63P, A71V, G73S, V82A and I84V Single Mutants.

A great challenge toward Acquired Immunodeficiency Syndrome (AIDS) treatment is to combat the HIV-1 virus. The major problem of drug resistance has kept the virus one step ahead of the medical community, and the call for more effective drugs remains as urgent as ever. Saquinavir, the first inhibitor against HIV-1 protease, offers the most extensive clinical data regarding resistance mutations. In this work, we examine L10I, G48V, L63P, A71V, G73S, V82A, and I84V single mutant HIV-1 PR strains in complexes with saquinavir to elucidate drug-protease interactions and dynamics. A comparative analysis of these mutations at the molecular level may lead to a deeper understanding of saquinavir resistance. The G48V mutation induces structural changes to the protease that reflect upon the drug's binding affinity, as shown by MM-PBSA and thermodynamic integration (TI) calculations (ΔΔGTI = 0.3 kcal/mol; ΔΔGMM-PBSA = 1.2 kcal/mol). It was shown that mutations, which increase the flexibility of the flaps (G48V, L63P, L10I) diminish binding. The preservation of hydrogen bonds of saquinavir with both the active site and flap residues in the wild-type and certain single mutants (A71V, V82A) is also crucial for effective inhibition. It was shown that mutations conferring major resistance (G48V, L63P, I84V) did not present these interactions. Finally, it was indicated that a water-mediated hydrogen bond between saquinavir and Asp29 in the active site (wild-type, A71V, G73S) facilitates a proper placement of the drug into the binding cavity that favors binding. Mutants lacking this interaction (G48V, V82A, I84V) demonstrated reduced binding affinities. This systematic and comparative study is a contribution to the elucidation of the drug resistance mechanism in HIV-1 PR.

Dual inhibitors for aspartic proteases HIV-1 PR and renin: advancements in AIDS-hypertension-diabetes linkage via molecular dynamics, inhibition assays, and binding free energy calculations.

Human immunodeficiency virus type 1 protease (HIV-1 PR) and renin are primary targets toward AIDS and hypertension therapies, respectively. Molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) free-energy calculations and inhibition assays for canagliflozin, an antidiabetic agent verified its effective binding to both proteins (ΔG(pred) = -9.1 kcal mol(-1) for canagliflozin-renin; K(i,exp)= 628 nM for canagliflozin-HIV-1 PR). Moreover, drugs aliskiren (a renin inhibitor) and darunavir (an HIV-1 PR inhibitor) showed high affinity for HIV-1 PR (K(i,exp)= 76.5 nM) and renin (K(i,pred)= 261 nM), respectively. Importantly, a high correlation was observed between experimental and predicted binding energies (r(2) = 0.92). This study suggests that canagliflozin, aliskiren, and darunavir may induce profound effects toward dual HIV-1 PR and renin inhibition. Since patients on highly active antiretroviral therapy (HAART) have a high risk of developing hypertension and diabetes, aliskiren-based or canagliflozin-based drug design against HIV-1 PR may eliminate these side-effects and also facilitate AIDS therapy.

Binding of novel fullerene inhibitors to HIV-1 protease: insight through molecular dynamics and molecular mechanics Poisson-Boltzmann surface area calculations.

The objectives of this study include the design of a series of novel fullerene-based inhibitors for HIV-1 protease (HIV-1 PR), by employing two strategies that can also be applied to the design of inhibitors for any other target. Additionally, the interactions which contribute to the observed exceptionally high binding free energies were analyzed. In particular, we investigated: (1) hydrogen bonding (H-bond) interactions between specific fullerene derivatives and the protease, (2) the regions of HIV-1 PR that play a significant role in binding, (3) protease changes upon binding and (4) various contributions to the binding free energy, in order to identify the most significant of them. This study has been performed by employing a docking technique, two 3D-QSAR models, molecular dynamics (MD) simulations and the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method. Our computed binding free energies are in satisfactory agreement with the experimental results. The suitability of specific fullerene derivatives as drug candidates was further enhanced, after ADMET (absorption, distribution, metabolism, excretion and toxicity) properties have been estimated to be promising. The outcomes of this study revealed important protein-ligand interaction patterns that may lead towards the development of novel, potent HIV-1 PR inhibitors.

Conformational Properties and Energetic Analysis of Aliskiren in Solution and Receptor Site.

Aliskiren is the first orally active, direct renin inhibitor to be approved for the treatment of hypertension. Its structure elucidation and conformational analysis were explored using 1D and 2D NMR spectroscopy, as well as random search and molecular dynamics (MD) simulations. For the first time, MD calculations have also been performed for aliskiren at the receptor site, in order to reveal its molecular basis of action. It is suggested that aliskiren binds in an extended conformation and is involved in several stabilizing hydrogen bonding interactions with binding cavity (Asp32/255, Gly34) and other binding-cavity (Arg74, Ser76, Tyr14) residues. Of paramount importance is the finding of a loop consisting of residues around Ser76 that determines the entrapping of aliskiren into the active site of renin. The details of this mechanism will be the subject of a subsequent study. Additionally molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) free energy calculations for the aliskiren-renin complex provided insight into the binding mode of aliskiren by identifying van der Waals and nonpolar contribution to solvation as the main components of favorable binding interactions.